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Abstract 

This paper focuses on observer based fault reconstruction for a class of nonlinear uncertain 
systems with Lipschitz nonlinearities. An adaptive-gain Super- Twisting (STW) observer is 
developed for observing the system states, where the adaptive law compensates the uncertainty 
in parameters. The inherent equivalent output error injection feature of STW algorithm is then 
used to reconstruct the fault signal. The performance of the proposed observer is validated 
through a Hardware-In-Loop (HIL) simulator which consists of a commercial twin screw 
compressor and a real time Polymer Electrolyte Membrane fuel cell emulation system. The 
simulation results illustrate the feasibility and effectiveness of the proposed approach for 
application to fuel cell systems. 

I. INTRODUCTION 

Many physical systems cannot be operated safely without reliable Fault Detection and 
Isolation (FDI) schemes in place. Such systems include Fuel Cells, which use Hydrogen 
as fuel, and system failure or mechanical faults can lead to catastrophic consequences. 
FDI is usually achieved by generating residual signals, obtained from the difference 
between the actual system outputs and their estimated values calculated from dynamic 
models. Such approach usually involves two steps: the first step is to decouple the faults 
of interest from uncertainties and the second step is to generate residual signals and detect 
faults by decision logic. Several practical techniques for these steps have been proposed 
in contemporary literature, for example geometric approaches [?], H -optimization 
technique [?], [?], observer based approaches (e.g. adaptive observers [?], [?], high gain 
observers [?], unknown input observers [?], [?], [?]). However, in active fault tolerant 
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control (FTC) systems [?], not only the fault needs to be detected and isolated but 
also needs to be estimated so that its effect can be compensated by reconfiguring the 
controller [?], [?]. Hence, there is need for fault reconstruction schemes which estimate 
the fault's shape and magnitude. 

Sliding mode technique is known for its insensitivity to external disturbances, high 
accuracy and finite time convergence. Sliding mode observers have been widely used for 
fault reconstruction. [?] proposed a fault reconstruction approach based on equivalent 
output error injection. In this method, the resulting residual signal can approximate the 
actuator fault to any required accuracy. Based on the work of [?], [?] proposed a sensor 
fault reconstruction method for well-modeled linear systems by considering a linear 
matrix inequality (LMI) observer design. This approach is impractical, as there is no 
explicit consideration of disturbance or uncertainty. To overcome this, [?] proposed an 
FDI scheme for a class of linear systems with uncertainty, using LMI for minimizing the 
L 2 gain between the uncertainty and the fault reconstruction signal. It should be noted 
that, only linear systems are considered in the above works, and only few works have 
been reported for nonlinear systems. [?] proposed a sliding mode observer based fault 
estimation approach for a class of nonlinear systems with uncertainties. [?] proposed a 
precise fault reconstruction scheme, based on equivalent output error injection, for a class 
of nonlinear systems with uncertainty. A sufficient condition based on LMI is presented 
for the existence and stability of a robust sliding mode observers. The limitation is that 
requires a strong structural condition of the distribution associated with uncertainties. 
This structural constraint is relaxed in [?], where the fault distribution vector and the 
structure matrix of the uncertainty are allowed to be functions of the system output and 
input. However, in these papers involving uncertainty and fault, most of all require that 
the uncertainty and fault are bounded with known bounds. Although in the job of [?], 
it is not necessary to known the bounds on uncertainty a priori, but it still requires that 
the boundary of fault signal is known. 

From the application point of view, fuel cell systems have also been the subject of 
many FDI studies. [?] developed a hydrogen leak detection method without the use of 
relative humidity sensors. [?] designed a fault diagnosis methodology for PEMFCs where 
the residuals are generated from the differences between the PEM simulator included 
with a set of typical faults and non-faulty fuel cell model. [?] proposed a Linear Parameter 
Varying (LPV) observer based fault detection approach for PEMFC, the residuals are 
defined as the differences between the process measurements and the outputs of a LPV 
Luenberger observer. The above works have not considered the design of state estimation, 
parameter identification and directly fault reconstruction, simultaneously. Even more, the 
models used in [?], [?] are obtained through a Jacobian linearization of the PEMFC 
nonlinear dynamic model around the optimal operating point. It is well known that the 
fuel cell system exhibits highly nonlinear dynamics and a fuel cell's efficiency is highly 
dependent on the operating conditions, such as temperature, humidity, and air flow. A set 
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of auxiliary elements (valves, compressor, sensors, etc.) are needed to make the fuel cell 
work at the optimal operating point. The problems of state estimation and fault detection 
arise because of incomplete knowledge of the parameter and states of the system. 

In this paper, fault reconstruction is studied for a class of nonlinear uncertain systems 
with Lipschitz nonlinearities. This class of systems sufficiently defines the nonlinear 
dynamics of the Fuel Cell system. The uncertain parameter is estimated by a simple 
adaptive update law. Estimate of the uncertain parameter is then injected into a Sec- 
ond Order Sliding Mode (SOSM) observer based on Super Twisting algorithm, which 
maintains a sliding motion even in the presence of faults. The faults are reconstructed 
by analyzing the equivalent output error injection [?], [?]. We extend the result of [?] 
to a class of nonlinear uncertain systems with Lipschitz nonlinearities. The difficulty of 
the a-priori knowledge of fault signal bound, as encountered in [?], is overcome by the 
gain adaptation with respect to observation errors. The gains of Super Twisting (STW) 
algorithm are allowed to adapt based on the 'quality' of the obtained sliding mode, and 
the relation of the gains are chosen to make the algorithm hold homogeneity. Lyapunov 
analysis of the system with the adaptive observer demonstrates that the algorithm estab- 
lishes ideal SOSM, as compared to a similar work by [?], in which real Sliding Mode 
is established. 

The Fault reconstruction scheme is applied on a PEMFC system which belongs to a 
class of nonlinear uncertain systems with Lipschitz nonlinearities. The nonlinear model 
of a fuel cell system, presented in [?], is used. This model has been validated in a large 
operating range. The adaptive-gain STW sliding mode observer estimates the system 
states (oxygen and nitrogen partial pressures and compress speed), and an adaptive 
update law is employed to estimate the stack current, which is considered as an uncertain 
parameter. This eliminates the need of an extra current sensor. The fault considered in 
this study is a mechanical failure in the air circuit, which results in an abnormal air 
flow. The fault signal is reconstructed by analyzing the equivalent output error injection, 
which is obtained on-line from comparisons between the measurements from the sensor 
installed in the real system and the outputs of the observer system. In order to verify the 
effectiveness of the proposed approach of the application to a fuel cell system, Hardware- 
In-Loop (HIL) tests are conducted on a test bench which consists of a commercial twin 
screw compressor and a real time fuel cell emulation system [?], [?], [?]. The test bench, 
which is presented in this work, disposes a HIL simulator consists of a real time fuel 
cell system and a twin screw compressor. These test results show that the fault detection 
scheme successfully reconstructs the faults. 
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II. PROBLEM FORMULATION 
Consider the following nonlinear system, 

x = Ax + g(x, u) + <fi(y, u)9 + u(y, u)f{t), 
V = Cx, 

where x E W 1 is the system state vector, u(t) E R m is the control input which is 
assumed to be known, y E W is the output vector. g(x, u) E W 1 is Lipschitz continuous, 
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and u>(y,u) E 



r are assumed to be some smooth and bounded 
functions with p > q + r. The unknown parameter vector 9 E W is assumed to be 
constant and f{t) E W is smooth fault signal vector, which satisfies 
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where pi, p 2 are some positive constants that might be known or unknown. 

is an observable pair, and there exists a linear coordinate trans- 
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where LOi(y, u) E 

System is described by the following equations in the new coordinate system, 

z =TAT~ x z + Tg(T~ x z, u) + T<p{y, u)9 + Tu(y, u)f(t), 

y =ct~ x z, 

By reordering the state variables, system (HJ can be rewritten as 

y =A n y + A 2 iZ2 + gifa, V, u) + 0i (y,u)6 + ui (y,u)f(t), 
z 2 =A 22 z 2 + A 21 y + g 2 (z 2 , y, u) + <p 2 (y, u)9, 

y =zi, 



(4) 
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III. ADAPTIVE HIGH ORDER SLIDING MODE OBSERVER DESIGN 

Now, we consider the problem of an adaptive-gain STW observer for system (HJ, in 
which the uncertain parameters are estimated with the help of an adaptive law. Then, 
using the equivalent output injection of the observer, we will develop the precise fault 
reconstruction method. The basic assumption on the class of nonlinear systems under 
consideration is as follows: 



Assumption 2. There exists a nonsingular matrix T G 







rxq 



l pxp , such that 

0<jxr 
$2(0, u) 



(7) 



where $i(y,u) G W xq , & 2 (y,u) G W xr are both nonsingular matrices. 

Remark 1. The main limitation in the Assumption ([7|) 15 that the matrices [4>i(y, u) uii(y, w)] 
must be block-diagonalizable by elementary row transformations. For simplicity, the case 
p = q + r is considered, the same results can be obtained in the case p > q + r [?]. 



Let z y = Ty, where T is defined in Assumption ©. Then, System © can be described 



by, 



z y =TA U T 1 z y + TA 21 z 2 + Tg 1 (y, z 2 , u) + 
z 2 =A 22 z 2 + A 21 y + g 2 (z 2 , y, u) + <p 2 (y, u)9, 

y =t~ x Zv, 



®i(y,u) 
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--A llV + A 21 z 2 + W gi (y, z 2 , u) + $i(y, u)9, 
--A 12 y + A 22 z 2 + W g2 (y, z 2 , u) + $ 2 (y, u)f(t), 



z 2 =A 21 y + A 22 z 2 + g 2 (y, z 2 , u) + (p 2 (y, u)9, 

I T 
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The observer is represented by the following dynamical system 

z yi =A u y + A 21 z 2 + W gi (y, z 2 ,u) + $i(y,u)0 + /i(e yi ), 



z. 



--A 12 y + A 22 z 2 + W 92 (y, z 2 ,u) + /x(e 



z 2 =A 21 y + A 22 z 2 + g 2 (y, z 2 , u) + <f) 2 (y, u)0, 



(ID 



{ y=T- l [ 



z yi Z V2 



where the continuous output error injection fx(s) is given by the super-twisting algorithm 
[?]: 

= X(t)\s\^ sign(s) + ip(s) 



<f(s) = a(t)sign(s) 



(12) 



with the gains A (t), a(t) are functions of time and are explained later on in Theorem 
3.3. The observation errors are then defined as e yi = z yi — z yi , e V2 = z V2 — z y2 ,e 2 = 
z 2 — z 2 , 9 = 9 — 9. The estimate of 9, denoted by 9, is given by the following adaptive 
law: 

9 = -K(y,u) (A u y + A 21 z 2 + W gi {yJ 2 ,u) + ^ l {y,u)§-Zy 1 ) , (13) 
where K(y,u) G W. qxq is a matrix design parameter which will be determined later. 



Remark 2. It can be seen that the adaptive law 4731) depends upon z Vl . A real time 
robust exact differentiator proposed in [?] can be used to estimate the time derivative 
of z yx in finite time. The differentiator has the following form 



(14) 



= - Ao-^o \ z o ~ z vi 1 2 s W n ( z o ~ %J + zi, 
= - a L sign(z - z yi ), 

where z and Z\ are the real time estimations of z yi and z yi respectively. The parameters 
of the differentiator X = l,a = 1.1 are suggested in [?]. L is the only parameter 
needs to be tuned according to the condition \z yi \ < L . 

It follows from (flOl) and (ITTl [T3l that the error dynamical equation is 

e 2 = A 22 e 2 + g 2 (y,z 2 ,z 2 ,u) + (p 2 (y,u)9, 
9 = -K(y,u) (A 21 e 2 + W gi {y,z 2 ,z 2 ,u)+$i(y,u){ 
-/i(e yi ) + A 21 e 2 + + W 9l (y,z 2 ,z 2 ,u) 



-fi(e V2 ) + A 22 e 2 + $ 2 (y,u)f(t) + W g2 (y, z 2 , z 2 ,u) 



(15) 
(16) 

(17) 
(18) 



where g 2 (y, z 2 , z 2 ,u) = g 2 (y,z 2 ,u) - g 2 (y, z 2 ,u),W gi {y, z 2 , z 2 ,u) = W gi (y,z 2 ,u) - 
W gi (y,z 2 ,u) and 
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W g2 (y, z 2 , Z2, u) = W g2 (y, z 2 ,u) — W g2 (y, z 2 ,u). Some assumptions are imposed upon 
the system: 

Assumption 3. The known nonlinear terms g 2 (y,Z2,u), W gi (y, z 2 ,u) and W 92 (y, z 2 ,u) 
are Lipschitz continuous with respect to z 2 i.e. 

Hflbd/, 32, «) - 9t.(V, Z2, U) || < -f 2 \\z2 ~ Z 2 \\, 

\\W gi (y,z 2 ,u) - W gi (y,z 2 ,u)\\ < "y gi \\z 2 - z 2 \\, (19) 
\\W g2 (y,z 2 ,u) - W g2 (y,z 2 ,u)\\ < i g2 \\z 2 - z 2 \\. 

where j gi , ■y 92 and 72 are the known Lipschitez constants for W gi (y, z 2 , u), W 92 (y, z 2 , u) 
and g 2 (j/, z 2 , u) respectively [?]. 

Assumption 4. Suppose that the Huritwz matrix A 2 2 satisfies the following Riccati 
equation 

A T 22 P X + P 1 A 22 + 1IP1P1 + 2/ n _ p + el n - p = 0, (20) 

which has a symmetric positive-definite solution Pi for some e > [?]. 

Assumption 5. Suppose that the matrix design parameter K(y, u) satisfies the following 
equation 

K(y, it)$i(y, it) + #T(y, u)K T (y, u) - -y 2 gi K(y, u)K T (y, u) - el q = 0, (21) 
for some e > 0. 

Now, We will first consider the stability of the error systems (fT5l) and (TT6l) . 
Theorem 1. Consider the systems rfTTl U6\) satisfying the Assumptions (TJ] ElanJUJ. Then, 



(22) 



systems 4771) and rtioD are exponentially stable if the matrix 

eln-p P\4>i{y, u) - A 21 K T (y, u) 



Q 



<f) 1 2 (y,u)P 1 -K(y,u)A2i el, 



'/ 



15 positive definite, where P\ and K(y, u) satisfies rf20l) and ([27]) respectively. 
Proof: A candidate Lyapunov function is chosen as 

V^J) = e 2 P P 1 e2 + 9 T 6, (23) 
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and the time derivative of V\ along the solution of the system (fT5l [T6l) is given by 

Vi =el{A T 22 P l + P 1 A 22 )e 2 + 2e^P 1 g 2 (v, *2, z 2 , u) + 2e T 2 P l <p 2 {y, u)9 - 2e T 2 A T 21 K T {y , u)~9 

- 29 T K(y, u)W gi (y, z 2 , z 2 , u) - 9 T (K(y, u)$x(y, u) + ®J(y, u)K T (y, u)) 9 
<e T 2 {A T 22 P x + P,A 22 + ^PxPi + 2I n . p )e 2 + 2e\ (PiMv, «) - A 21 K T (y, u)) 9 

- 9 T (K(y, u)^(y, u) + $J(y, u)K T (y, u) - j 2 gi K(y, u)K T (y, u)) 9 



-[4 9 T ]Q 



(24) 



Hence, It follows that e 2 , 9 will converge to zero exponentially. ■ 
Before proving the convergence of the systems (fTTt and (fTST ). we will show some 
conditions. Theorem [T] shows that lim e 2 (t) = and lim 9{t) = 0. Consequently, the 

t— >oo t— >oo 

errors e 2 ,9 and its derivatives e 2 , 9 are bounded. 

Assumption 6. The time derivative of functions $i(y,u), § 2 (y,u), W 9l (y,z, u) and 
W g2 (y, z, u) are bounded. 

Therefore, under the Assumption |6] the time derivative of the nonlinear terms in the 
error dynamics ([171 [T8T > are bounded. That is, 



j t (A 22 e 2 + $ 2 (//. //)/(/) - \V,J/j. y>. i- 2 . if . 
j [A 21 e 2 + ^ 1 {y,u)9 + W gi {y,z 2 ,z 2 ,u] 



< 



Xi, 



(25) 



< 



where xi an d X2 are some unknown positive constants. 

In what follows, the objective is to prove the finite time convergence of the systems 
(fTTb and CLE) based on adaptive-gain STW algorithms [?], [?]. 

Theorem 2. Consider the system ( 1771 ITffl) with the following adaptation laws, 

\{t) = 2y/L(jj, 

a(t) = 4L(t), 

where the dynamic of the positive time varying function L(t) is given by 

7^0 



(26) 



L{t) 



if 
else 



(27) 



where k > is a positive design constant. Then, all trajectories of the system ( 1771 U~8\) 
converge to zero in finite time. 
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The proof of this theorem is similar to [?]. 



Remark 3. The advantage of the law A26v is that the uncertainty affecting the system 
can be attenuated by —, where L is a positive constant. Consider the structure of the 
system 4771) and rf78l) . 

de , . ,i . , . 

— = - \ \e\2sign(e) + <p, 

f (28) 

dip . , v . . 

— = - a sign(e) + u, 

where v denotes the uncertainty which affects the system 4771) and ( EH ). 
Then according to (l26l) . the time-scale is modified as 

dr = Ldt, e = Le. (29) 

then, for every fixed value of L, the system ft28\) is equivalent to the following system: 

de . . ... 

— = - \ \e\2sign(e) + <p, 

f . (30) 

dip v 

— = - a sign(e) + -, 
dr L 

Now, it follows from a comparison between the systems rf2g|) and 001) . that the uncertainty 
effects on the system (TJOj) will be attenuated by — for every fixed positive value of L with 

respect to the new time-scale r, without changing the system structure. In our paper, 
A = 2 and a = 4 are chosen. 

Theorems \T\ and |2] have shown that systems (fTTT) and (fi~3T) are an asymptotic state 
observer and uncertain parameter observer for the system (llOl) respectively. In the next 
Section, we will develop the fault reconstruction approach based on those two observers. 

IV. Fault Reconstruction 

The fault signal f(t) will be reconstructed based on the proposed observer by using 
an equivalent output error injection which can be obtained as soon as the sliding surface 
is reached. 

It follows from Theorem |2] that e yi , e V2 , e yi , e V2 in (flTT) and ([TBI are driven to zero in 
finite time. Then, the equivalent output error injection are obtained directly 

H{e y2 ) = A 22 e 2 + $> 2 (y,u)f(t) + W g2 (y,z 2 ,z 2 ,u). (31) 

From Assumption |2l & 2 (y,u) is a nonsingular matrix, then the estimation of f(t) can 
be constructed as 

f(t) = <S> 2 1 (y,u)n(e y2 ). (32) 
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Theorem 3. If the conditions of Theorems ([7J [2]) are satisfied, then f(t) defined in rfJ2 
is a reconstruction of the fault fit) since 



lim 

t— >oo 



0. 



/(*) - fit) 

Proof: It follows from CD} and ([32]) that 

|/(t) - /(t)|| = \\^\y^) (^2262 + W 92 

< \\<5> 2 \y,u)A 22 \\ ||e 2 || +7 32 ||$ 2 ^y, u) || ||e 2 
It follows from Theorem CD that 



lim 

t— >oo 



fit) - fit) 



0. 



(33) 



(34) 



(35) 



Hence, Theorem |3] is proven. ■ 

V. Application to PEM Fuel Cell systems 

In this section, we demonstrate the observer design process using the reduced-order 
model of a fuel cell system which has been verified in [?] experimentally. The dynamic 
model of the fuel cell system is given as, 



where the vector F(x) 
h{x)=- 

/s(*) 
Uix) 



x =F(x) + g u u + g£ + g f f, 

y =h(x), 

[fiix) hix) h{x) U{x)] T is 

cs{xi - c 2 )^ix x ) 



(36) 



(ci + c 8 )(xi - x 4 ) 

c 3 x 2 tp(x 1 



KXi 



- c 8 {x 1 - 

- C9X3 - 



• x^ 

ClO 

£3 



, X4 ■ 



C12 



- 1 



jl + Cis 









^3(^3), 

x [^3(^3) + ci 6 (xi - 2 4 )] 



(37) 



and the input vectors are g u = [0 C13 0] T , ^= [— C4 0] J , = [0 C5]' 
The system states 2 = [xi, x 2 , £3, x 4 ] T , where sci and s 2 are the total pressure at 
the cathode and nitrogen partial pressures respectively, x 3 is the angular speed of the 
compressure, and x 4 is the air pressure in the supply manifold; The measurements of 
the system are y = h(x) = [yi y 2 ] = [xi xA . The stack current £ is considered as 
an uncertain parameter 9 in our observer design, which eliminates the need of adding 
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current sensors in the system. The control input u represents the motor's quadratic current 
component is controlled using a real time static feed-forward controller. 

A failure of the fuel cell air circuit is considered, i.e. the pipe connecting the air- 
feed compressor to the fuel cell cathode explodes suddenly. This can be detected as a 
variation in supply manifold pressure dynamics [?]. The fault signal f(t) appears in the 
output channel z y2 , 



f(t) 



3 x 10 3 kg/sec, if t > 50 sec 
0, else 



(38) 



The function ip(xi) is the total flow rate at the cathode exit 

V>(xi) : 



CyjX\ 



Cll 



Xi 



Cll 
X\ 



C20X\, 



if > C19 

X\ 

., C11 , 

IJ < C19 

Xi 



(39) 



and h(x 3 ) is the mass flow rate of a twin screw compressor, where h(x 3 ) = Ax?,. All the 
parameters q, kf are positive and depend on the physical values of the fuel cell [?](See 
Appendix A). 

In order to design the proposed observer for the fuel cell system. Let define z yi := x\, 
z y2 := £4, z 2 := [x2 £3] and 9 := £. Then, system (|36l ) is described as the form of 



Z-2 



(c 1+ c g )fe,-z w )- C3fa '- C ^ fa - ) -cA 



KZ 



Ul 



jl + Cis 









x [h 3 (D 2 z 2 ) + c 16 (z yi - z y2 ) - f(t)] 



-H " 


Z2 + 








z Vi 






Zy 2 ^ 







-c 8 c 8 


A21 



y + 



ClO 



D 2 z 2 



D lZ2 \ H 

z, XC12 



HZ 



V2 
Cll 



HI 



h 3 {D 2 z 2 ) + c 13 u 



92(y,z2,u) 



where D 1 = [l 0] , D 2 = [0 l] , W gi (y, z 2 , 



u 



c 3 {z 



C2)i>{Z; 



(40) 



VI 1 



KZ 



W g2 (y,z 2 ,u) := C14C15 



^V2 

Cll 



Hi 



(h 3 (D 2 z 2 ) + c 16 (z yi - ZyJ) , $i(y, u) := -c 4 , $2(2/, «) - 
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C5, 4>2(y, u) — and the design gain parameter H is chosen to satisfy the Riccati equation 
(l20l) . The fault signal is weighted by c 5 , modeled as 

C12 



^5 



-Ci4 



1 + 



(41) 



The adaptive-gain STW observer for the system (l40l is designed as the form (TTTT IT3T) 



-:vi 



'1)2 



Z2 



{ | \ / \ C 3(^i — c 2 )i>{z yi ) a s ( \ 
(d + cs)(^ - 2^, ) * 2 C4 + //(e ai ), 



K,2 



Cl4 S 1 + Cl 5 



'^2 \C12 

Cii 



-c 9 D 2 z 2 



-c 8 (z 

ClO 



J/1 " ^2/2/ 



x [/i 3 ( J D 2 ^) + c 16 (z yi - z y2 )] + //(e^), 
c 3 D 1 z 2 i'(z yi ) 



HZ 



D 2 z 2 



fZ V2 \C12 _ I 



III 



h 3 {D 2 z 2 ) + c 13 u 



(42) 



'in 



-))■> 



and 



6 = -K I cj + (d + c 8 )(z 



ui 



Zy 2 ) -t -t Z yi 



HZ 



!)1 



(43) 



the adaptive-gains of the STW algorithm n(e yi ), [i(e y2 ) are designed according to (1261 
. The value of z yi is obtained from the robust exact finite time differentiator (TBI) in 



[?]. From (l32l . the fault signal /(£) is estimated as / 



n(e. 



1)2 I 



Remark 4. The Assumptions (\3\6\) are satisfied by the functions W gi (y, z 2 , u), W g2 (y, z 2 , u), 
$x(y,u), & 2 (y,u) and g 2 (y,z 2 ,u). The Riccati equation in the Assumption \4\ will be 
satisfied by appropriate value of the design gain H > 0. The Assumption \5\ is also 
satisfied for some e > 0, since the equation ([27]) is simplified into a scalar equation. 



VI. Experimental Results 

Experiments have been performed on a Hardware-In-Loop (HIL) test bench shown 
in Fig. 1, which consists of a real time emulated fuel cell system and a twin screw 
compressor. This emulated PEMFC stack is a 33 kW fuel cell composed of 90 cells in 
series, which provides the cathode pressure and supply manifold pressure as outputs. 
The twin screw compressor consists of two helical rotors which are coupled directly 
to its motor. Air intake is at the opposite side of the mechanical transmission and the 
output pressure is regulated by a servo valve. It has a flow rate margin 0-0.1 kg/s at 
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a maximum velocity of 12000 RPM, and is driven by a permanent synchronous motor 
(PMSM). The lubrication system of this compressor is specifically designed to prevent 
contamination of the air from lubricating oil. 

The test bench is controlled by the National Instruments CompactRIO realtime con- 
troller and data acquisition system. The motor quadratic current calculated by the real 
time controller in (d, q) coordinate is transformed by an inverter between the Com- 
pactRIO and the compressor to 3-phase coordinate in order to control the PMSM. The 
fault reconstruction scheme and HIL simulation structure is shown in Fig. 2. The nominal 
value of the parameters for the simulation are shown in Table 1 . 

The air mass-flow measured at the output of the compressor is fed to the real time 
fuel cell emulation system, programmed in Labview. The stack current shown in Fig. 7 
is varied between 100A and 450A in order to deal with different load variation which 
corresponds the flow rate variation between Og/s and 28g/s. The detail of the real time 
controller is available in [?]. 

Figs. 3-6 show the states estimation of the system (T3~6l ). Fig. 7 shows that the adaptive 
law gives a good estimate for the stack current which is considered as an unknown 
parameter 6. Based on the estimates of the proposed observer, the fault signal is re- 
constructed faithfully as shown in Fig. 8. The time history of the adaptive-gain L(t) is 
shown in Fig. 9, where the convergence of the proposed observer is ensured. 

VII. CONCLUSION 

This paper has proposed a robust fault reconstruction method for a class of nonlinear 
uncertain systems with Lipschitz nonlinearities based on an adaptive STW sliding mode 
observer. An adaptive update law has been given to identify the uncertain parameters. 
An adaptive-gain STW observer is proposed to estimate the system state variables expo- 
nentially even in the presence of uncertain parameters and fault signals without requiring 
any information on the boundaries of the fault and its time derivative. Furthermore, the 
obtained equivalent output error injection was then used to reconstruct the possible faults 
in the system. The proposed fault reconstruction approach was successfully implemented 
on a Hardware-In-Loop test bench which consists of a commercial twin screw compressor 
and a real time fuel cell emulation system. The fault signal represents as a failure of 
the fuel cell air circuit, i.e. the pipe connecting the air-feed compressor to the fuel cell 
cathode explodes suddenly, was reconstructed precisely. The experimental results have 
shown that the proposed approach is effective and feasible. In future, other faults that 
affect fuel cell performance, such as drying or flooding at the cell stack and starvation 
will also be included. 
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Twin screw 




Fig. 1. Test Bench 




Fig. 2. 



Fault Reconstruction Scheme implemented on a Hardware-In-Loop Simulator 
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TABLE I 

Parameters used in the simulation 



Symbol 


Parameter 


Value 


n 


Number of cells in fuel cell stack 


90 


h 


Universal gas constant 


0.314 J/(mol Kj 




Air gas constant 


zoo.y j/^Kg ps.) 


patm 


Atmospheric pressure 


I.UIjZj X1U ra 


T 

aim 


Atmospheric temperature 


ZVo.lJ JS. 


J-fc 


Temperature of the fuel cell 


333. is. 


1 


faraday constant 


704OJ wmoi 


M a 


Air molar mass 


28.9o44 x 1U kg/mol 


Mo 2 


Oxygen molar mass 


TO , , 1 A — 3 1 / 1 

32x10 kg/mol 


M N2 


Nitrogen molar mass 


28xl0~ 3 kg/mol 


M v 


Vapor molar mass 


18.02x10 kg/mol 


<~-D 


Discharge of the nozzle 


U.UU3o 


A 

At 


Operating area of the nozzle 


0.00138 m 


7 


Ratio of specific heats of air 


1 i 
1 .4 




1 nrnt^fp cent" 1 r\&rtt q 
\^.UllipiCaSUI lilCI Lid. 


671 Q yin -5 ktr m 2 

0/1.7 AlU Kg / / 1 


f 
J 


IVIotor friction 


u.uuiju v/^rati/s^ 


u 
lit 


Motor constant 


U.31 IN 111/ A 




Constant pressure specific heat of air 




l/cp 


rnmnrpisnr pffinpnrv 


80% 




Motor mechanical efficiency 


98% 


v ca 


Cathode volume 


0.0015 m 3 




Supply manifold volume 


0.003 m 3 


Vcpr / tr 


Compressor volume per turn 


5xl0~ 4 m 3 //r 




Cathode inlet orifice constant 


0.3629xlO~ 5 fcp/(Pas) 


fcca,out 


Cathode outlet orifice constant 


0.76xlO~ 4 fcg/(Pas) 


Pa 


Air density 


1.23 kg/m 3 


,ea,in 


Oxygen mass fraction 


0.23 




Smoothing filter time constant 


0.01 s 


L(0) 


Initial value of the adaptive gain L(t) 


5000 


k 


Design parameter in i27\ 


500 



0.15 - 



c 



10 20 30 40 50 60 70 

Time (s) 



90 100 110 



Fig. 3. Oxygen partial pressure estimation 
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Fig. 4. Nitrogen partial pressure estimation 
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Fig. 5. Compressure speed estimation 
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Fig. 6. Supply manifold pressure estimation 
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Fig. 7. Adaptive estimation of the stack current £ 
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Fig. 8. Fault reconstruction / and / 
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Fig. 9. 



Adaptive law of L(t) 
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Appendix 

RT s tkca,in j Xo2,atm \ p 

M 02 V ca \l + u atm J 
RT st _ RT st n 

RT s tk ca ^i n /l Xc>2,atm 

C8 ■ 



M N2 V ca V l + Uatm , 

c 9 — — j — ; c io — - j 

t> _ 7- 1 

Cn — r atm , C12 — 

7_ 

Vcmh 'jRTatm 

C13 = — ; — ; Cu 



Cl5 5 Cxg ^ca,ira 

cdAt n^r i 

Cl7 = W r! Gig = - 

\/RJ\ t S 7-1 7 
1 / 2 \^ 

^4 = —Tlv-cVcpr ItrPa, c 19 = I ~~T 

27T \7 + 1 

C D A t i / 2 \ ^ 



c 20 = ^ — 7 2 < . n 



